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Abstract 

We discuss the result of a recent quenched lattice calculation of the K — ► it vector form factor 
at zero-momentum transfer, relevant for the determination of \V US \ from K — * niv decays. Using 
suitable double ratios of three-point correlation functions, we show that it is possible to calculate 
this quantity at the percent-level precision. The leading quenched effects are corrected for by means 
of quenched chiral perturbation theory. The final result, f+ n (0) = 0.960 ± 0.005 sta t ± 0.007 ayst , 
turns out to be in good agreement with the old quark model estimate made by Leutwyler and 
Roos. In this paper, we discuss the phenomenological impact of the lattice result for the extraction 
of \V US \, by updating the analysis of K — ► iriv decays with the most recent experimental data. We 
also present a preliminary lattice study of hyperon S~ — > nlv decays, based on a similar strategy. 



1 Based on talks given at: DA&NE 2004: Physics at meson factories, Laboratori Nazionali di Frascati 
(Italy), June 7-11, 2004; VIII International Conference on "Electron-Nucleus Scattering", Marciana Marina 
(Italy), June 21-25, 2004; Lattice 2004, Fermi National Accelerator Laboratory Batavia, Illinois (USA), June 
21-26, 2004; ICHEP 2004, Beijing (China), August 16-22, 2004. 



1 Introduction 



Semileptonic decays of kaons and hyperons are of great phenomenological interest, as they 
provide a determination of the CKM matrix element |V^ S | ll.J. Whereas the most precise 
estimate of \V US \ is presently obtained from K — > niv (K^) decays, the studies of hyperon 
decays provide an important, independent approach [2|. Other methods to extract |V^ S | are 
based on leptonic kaon decays [3J H] and r decays [5] . 

The analysis of the experimental data on K — > nlv 6J and hyperon |5j decays 

gives access to the quantity \V US \ ■ /y(0), where /y(0) is the vector form factor at zero- 
momentum transfer for each decay, i.e. /y(0) = /+(0) for K i3 decays and /y(0) = /i(0) 
for hyperon decays. Vector current conservation guarantees that, in the SU(3)-symmetric 
limit, /+(0) = 1 and /i(0) = —1 for the S~ — > niv transition. A good theoretical control on 
these decays is obtained via the Ademollo-Gatto (AG) theorem |7j, which states that /y(0) 
is renormalized only by terms of at least second order in the breaking of the SU(3)-flavor 
symmetry. The estimate of the difference of /y(0) from its SU(3)-symmetric value represents 
the main source of theoretical error and, in the case of decays, it dominates the overall 
uncertainty in the determination of \V US \. 

The amount of SU(3) breaking due to light quark masses can be investigated in the frame- 
work of Chiral Perturbation Theory (ChPT), where the form factors can be systematically 
expanded as f v (0) = 1 + f 2 + fs + ■ ■ with /„ = C(M"/A") and A x ~ 1 GeV. In the meson 
sector, only the even terms of this expansion appear. Moreover, thanks to the AG theorem, 
the first non-trivial term can be computed unambiguously in terms of Mk, M n and f n : 
f'2 = —0.023 in the K° — > 7r~ case jH]. At higher orders, the estimate becomes more difficult 
due to the presence of local effective operators whose coupling are essentially unknown. The 
next-to-leading order correction, / 4 , was evaluated many years ago by Leutwyler and Roos 
[Bj. By employing a general parameterization of the SU(3) breaking effects in the meson 
wave functions, they obtained / 4 = —0.016 ± 0.008. This result still represents the reference 
value adopted by the PDG [6 . It yields the estimate /+(0) = 0.961 ± 0.008. 

The complete ChPT calculation of / 4 for K — > 71 decays has been recently performed 
Elj. The whole result is the sum of a loop amplitude (including both genuine two- 
loop terms and one-loop diagrams with 0(p A ) couplings) and a local term that involves 
a combination of the (unknown) 0(p 6 ) chiral coefficients. The loop amplitude exhibits a 
large scale dependence ^T], which is compensated for by the corresponding variation of the 
0(p 6 ) couplings, to yield a scale-independent / 4 . An important observation by Bijnens and 
Talavera [Sj is that the combination of 0(p 6 ) constants appearing in / 4 could be constrained 
in principle by experimental data on the slope and curvature of the scalar form factor. The 
required level of experimental precision, however, is far from what can be presently achieved. 
Thus, for the time being, the only analytic estimates of the form factor are provided by either 
the Leutwyler- Roos result or other model dependent calculations [HHI2], and the large scale 
dependence of the 0(p 6 ) loop amplitude [TT] seems to indicate that the ±0.008 error might 
be underestimated. 

As for hyperon decays, the chiral expansion converges less rapidly (also the odd f\ are 
permitted) and the theoretical estimate is quite involved, already at the lowest order. The 
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corrections f2 can be unambiguously determined only in the limit in which baryons are 
treated as heavy particles and the octet-decuplet splitting is treated as a large scale [T3]. 
This splitting, however, is just of the order of the .fT-meson mass, and it was shown that the 
intermediate decuplet states can partially cancel the (formally leading) octet contribution 
|14j . Nothing is known about the higher orders terms. In addition, relativistic corrections, 
proportional to inverse powers of the "heavy-baryon" mass m#, have not been investigated 
so far. Since numerically mg ~ 47r/„., these contributions are expected to be of the same 
order of magnitude of the chiral corrections. Therefore, SU(3)-breaking effects in /i(0) are 
mainly estimated at present by using phenomenological models ^Hj, which yield corrections 
of the order of a few percent. 

It is clear that a reliable estimate of \ V US \, from either kaon or hyperon semileptonic decays, 
can be only based on a first principle determination of the SU(3)-breaking effects. Such a 
determination can be provided by lattice QCD, a non-perturbative approach based only on 
the fundamental theory. The first lattice calculation, with significant accuracy, of the SU(3)- 
breaking effects in K — > it transitions has been recently performed PH| (see also [T7]-|2(j|). 
A new strategy has been proposed and successfully applied in the quenched approximation 
in order to reach the challenging goal of about 1% error on /+(0). 

The aim of this paper is twofold. We first present the quenched lattice result for 
decays ^B] and discuss its impact on the extraction of |V^ S |, particularly in the light of the 
most recent experimental determinations. Second, we present the results of an exploratory 
study of hyperon decays based on the same strategy |2T] . This study has shown that the 
same level of statistical accuracy reached for decays can be also achieved for hyperon 
decays. In this latter case, however, further investigations within ChPT are required in 
order to better control the extrapolation to the physical meson masses, necessary in lattice 
calculations, and to reach a level of systematic accuracy comparable to the one achieved in 
the Ki 3 case. 

2 Semileptonic K —> tt£u decay 

2.1 Lattice result and phenomenological implications 

The K° ^ 7i~ form factors of the weak vector current = sj^u are defined as 

(7r(p')\V,\K(p)) = / + (g 2 )(p + p')M + f-tf)(p-P% , (1) 
where q = p — p' . As usual, f-(q 2 ) can be expressed in terms of the scalar form factor 

/ (g 2 ) = f + (q 2 ) + M /_ M J -(q 2 ) , (2) 

so that by construction /o(0) = /+(0). 

The procedure developed in Ref. ^Bj to reach the challenging goal of a about 1% error 
on /+(0), is based on the following three main steps: 

1) precise evaluation of the scalar form factor fo(q 2 ) at q 2 = g 2 ^ = (M K — M w ) 2 ; 
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o u yz j \ / 


6™ 

3-body 


(%) 

full 




2.31 ± 0.22 


-0.35 ± 0.16 


-0.10 ± 0.16 







+0.30 ± 0.10 


+0.55 ± 0.10 


K + 


2.31 ± 0.22 


-0.05 ± 0.20 


+0.20 ± 0.20 


K° 





+0.55 ± 0.20 


+0.80 ± 0.20 



Table 1: Summary of the isospin-breaking factors entering Eq. Q): 5^ [3 body] denotes the 
correction for the inclusive rate excluding the radiative events outside the Dalitz Plot 
Fll\ IHf; 5^ [full] denotes the correction for the fully inclusive rate ffltf . 

2) extrapolation of /o(<?Lx) to /o(0) = /+(0); 

3) subtraction of leading chiral logs and extrapolation of /+(0) to the physical meson 
masses. 

These steps are described in the next subsections and details of the simulation can be found 
in Ref. [TSj. Here we present our final result for the form factor at zero momentum transfer 
and discuss its phenomenological implications, particularly in the light of the most recent 
experimental results. 

By following the procedure outlined above, we obtain the result [16J 

ff*~ (0) = 0.960 ± 0.005 stat ± 0.007 syst , (3) 

where the systematic error does not include the estimate of quenched effects beyond 0(p 4 ). 
Eq. © compares well with the value ff w "(0) = 0.961 ± 0.008 quoted by the PDG ;6J and 
based on the Leutwyler-Roos quark model estimate of f± [Hj, thus putting the evaluation of 
the K — > 7r vector form factor on a firmer theoretical basis. 

The result in (|5J) . as well as the Leutwyler-Roos estimate, refers to the K° — > n~ form 
factor in absence of electromagnetic corrections. This quantity is a convenient and well- 
defined common normalization for all the four physical photon-inclusive decay widths 
(K = {K + , K }, I = {e, //}). By using this definition we can write 

r(K a ) = ^^iv^s^ff^m 2 ^^) [i + ^ (2) + & , (4) 

where Gf = 1.1664 x 10~ 5 GeV -2 , Ck = 1 (2 -1 / 2 ) for neutral (charged) kaon decays, 
Ik(\%) is the kinematical integral evaluated in the absence of electromagnetic corrections 
and S cw = 1.0232 is the universal short- distance electromagnetic correction renormalized at 
the scale \x = M p . The residual long-distance component of the electromagnetic corrections 
(which also need to be evaluated at \i = M p in order to cancel the scale dependence of 5 CW ) 
is encoded in 5fjj£. Note that 5£* depends both on the channel and on possible experimental 
cuts applied to the soft-photon radiation. Finally, the strong isospin-breaking correction due 
to m u ^ m d is encoded in 5f^, 2 N. The values of these channel-dependent corrections, which 
can be extracted by the recent analysis of Refs. [TT| |22l I2H] , are reported in Table HJ 
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Inverting Eq. (0J) we get 
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which could lead to four independent determination of |T4s| • /+ 71 (0) from the four 
modes. Three of these determinations have recently been updated by several experiments. 
A collection of these recent results is presented in Table El and illustrated in Fig. ^ In 
this plot we also show the averages of the old results quoted by the PDG. The new data 
consistently point toward a value of \V US \ ■ f^° n (0) substantially larger with respect to the 
one derived by the old PDG results. The value of \V US \ extracted from these new data, 
combined with our estimate of f+ 77 (0), is 

\V US \ = 0.2250 ± 0.0020 /+(0) ± 0.0005 cxp = 0.2250 ± 0.0021 . (6) 

Interestingly, this result is perfectly consistent with the value derived by CKM unitarity, 
thus solving a long-standing problem (see e.g. Ref. [30]). Using the accurate determination 
of \V u d\ from nuclear + — > + and nucleon beta decays, \V U d\ = 0.9740 ± 0.0005 [2H], one 
finds indeed |K s | unit ' - v 71 - lKd| 2 = 0.2265 ± 0.0022. 

Given our close agreement with the Leutwyler-Roos result, the solution of the CKM 
unitarity problem perfectly holds also if their estimate of f+ v (0) is adopted. On the other 
hand, the consistency is substantially worse with the value f K w (0) = 0.976 ± 0.010 of 
Ref. p. It is also worth to mention that the estimate f K ° w ~(0) = 0.981 ± 0.010 of Ref. [H] 
is not independent from the one of Ref. 0, and should not be used in conjunction with 
Eq. (J3J) and/or compared with our result. The difference between Ref. [TT] and most other 
analyses is due to a partial inclusion of the electromagnetic corrections in the form factor 
(concerning SU(3) breaking effects, they simply use the results of Ref. In our treatment, 
as well as in the analysis performed by the KTeV Collaboration all the electromagnetic 
corrections are explicitly factorized out. This approach has the advantage of separating the 
hard theoretical problem of estimating SU(3) breaking effects (identical for the four decay 
modes), from the less severe (but channel-dependent) issue of estimating electromagnetic 
corrections. 



2.2 Evaluation of f (q 2 ) at q 2 = q max = (M K - M^) 2 

In the next subsections we discuss the strategy of the lattice calculation of the vector form 
factor. A more exhaustive discussion and all details of the numerical simulation can be found 
in Ref. fTo]. 

By following a procedure originally proposed in Ref. [3T] to study heavy-light form factors, 
the scalar form factor fo(q 2 ) can be calculated very efficiently at q 2 = g^ ax = (M K — M n ) 2 
(i.e. p = p' = q = 0) from the double ratio of three-point correlation functions with both 
mesons at rest: 

ftU £ \ _ C^ttx, ty, 0, 0) Cg K (t x , ty, 0, 0) 
(t X , ty, 0, 0) C^ n (t X , ty, 0, 0) 
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5 13(10) 


5.14(6) 


27 01(9) 


27 02(25) 


\v us \ ■ f« u «- (0) 


0.2190(23) 


0.2192(15) 


0.2157(11) 

21 e3 

KLOE^W 


0.2157(14) 




KTeV[K L ] 


NA48[X L ] 


KL0E[iv~ 5 ]W 


BR(%) 


40.67(11) 


40.10(45) 


39.85(35) 


0.0709(11) 


\v us \-f^-(o) 


0.2160(9) 


0.2145(15) 


0.2138(13) 


0.2164(17) 



Table 2: Summary of recent results on decays obtained by BNL-E865 KTeV 
NA48 Wty and KLOE Jlffi and corresponding values of \V US \ ■ /+(0). The results marked by 
(*) are still preliminary. The values of \V us \f+ 77 (0) have been computed using the PDG 
values for the lifetimes [6l and using, in all channels, the slopes fitted by KTe V with a polar 
fit 



0.225 



0.222 



0.219 



0.216 



0.213 



0.210 



0.207 



BNL 



NA48 



LATTICE + UN1TARITY 



1 



t 



KTeV 



KLOE^'-Kg 



KLOE 



KLOE 



KTeV 



EXP^ 



6 



K 



e3 



e3 l£ M-3 



Figure 1: Experimental results for \ V US \ ■ /+(0). The gray (darker) band indicates the average 
of the new experimental results, \V US \ ■ /+(0) = 0.2160 ± 0.0005, whereas the yellow (lighter) 
band represents the unitarity prediction combined with our determination of the vector form 
factor: \V us \ mit ■ / + (0) = 0.2175 ± 0.0029 (see text). The lower (black) points without labels 
are the old PDG values. 
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where the three-point correlation function for the K — > 7r transition is defined as 



with Vfj, being the renormalized lattice weak vector current and 0\ = J75U, 0^ K = d^s the 
meson interpolating fields. 

When the vector current and the two interpolating fields are separated far enough from 
each other, the contribution of the ground states dominates the correlation functions, yielding 

p /, , x (n\sj u\K) (K\u-f s\n) 2 2 (M K + M n ) 2 
K {t xj ty) > - [/olg max jJ TT7 77 • (9 J 



— * 00 
tx) 



Z^^s^K) {h\u 1o u\tt) wu ~' j 4M^M ot 



Thus from this ratio the scalar form factor at maximum q 2 can be directly determined. 

There are several crucial advantages in the use of the double ratio (J7J). First, there is a 
large reduction of statistical uncertainties, because fluctuations in the numerator and the 
denominator are highly correlated. Second, the matrix elements of the meson sources cancel 
between numerator and denominator. Third, the double ratio is independent from both 
the renormalization constant and the 0(a)-improvement coefficient of the vector current, 
where a is the lattice spacing. Therefore Rq suffers from discretization effects only at order 
a 2 . Finally, the double ratio is equal to unity in the SU(3)-symmetric limit at all orders in 
a. Thus the deviation of Rq from unity depends on the physical SU(3) breaking effects on 
/o(?max) as wen as 011 discretization errors, but the latter are at least of order a 2 {m s — mi) 2 . 
In other words, our result for /o(<7max) * s n °t affected by the whole discretization error on the 
three-point correlation function, but only by its smaller SU(3)-breaking part. Similar con- 
siderations apply also to the quenching error, because the double ratio R (t x ,t y ) is correctly 
normalized to unity in the SU(3)-symmetric limit also in the quenched approximation. 

From the plateau of the double ratio Q the values of /b(<?max) are determined with a 
statistical uncertainty smaller than 0.1%, as it is illustrated in Fig. 



2.3 Extrapolation of f (qLx) to f o(0) = f+(0) 

The extrapolation of /o(<? 2 ) from g^ ax to q 2 = requires the knowledge of the slope of /o, 
and thus the study of the g 2 -dependence of the scalar form factor. An important remark is 
that, in order to get /o(0) at the percent level, the precision required for the slope can be 
much lower (~ 30%), because the values of g^ ax in the numerical simulation can be chosen 
very close to q 2 = 0. 

For each set of quark masses two- and three-point correlation functions have been calcu- 
lated for mesons with various momenta in order to extract the q 2 dependence of both fo(q 2 ) 
and f + (q 2 ). The latter turns out to be well determined on the lattice with a statistical 
error of about 5 20% and its g 2 -dependence is very well described by a pole-dominance fit, 
f+{q 2 ) = /+(0)/[l — A + q 2 ]. On the contrary, for the scalar form factor the uncertainties are 
about 5 times larger. As explained in Ref. ^Hl the precision in the extraction of fo{q 2 ) can 
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Figure 2: Values of fo(q 2 VAX ), obtained in Ref. fTb^ using Eq. (OJ), versus the SU (3) -breaking 
parameter a 2 AM 2 = a 2 [M 2 K - M 2 ). 



be significantly improved by constructing new suitable double ratios 

P u + \ _ Cf^itx, t y ,p,p') CQ K (t x , t y ,p,p') 
Ax,y) ~ C«*(t x ,t y ,p,p') C« K (t x ,t y ,p,p>) ' liUj 

(with i = 1,2,3) from which a determination of the ratio of the form factors fo{q 2 ) / f + (q 2 ) 
is obtained. The advantages of the double ratios (ITUjl are similar to those already pointed 
out for the double ratio (J7J), namely: i) a large reduction of statistical fluctuations; ii) 
the independence of the improved renormalization constant of the vector current, and iii) 
Ri — > 1 in the SU(3)-symmetric limit. We stress that the introduction of the matrix elements 
of degenerate mesons in Eq. (jTUj) is crucial to largely reduce statistical fluctuations, because 
it compensates the different fluctuations of the matrix elements of the spatial and time 
components of the weak vector current. 

Thanks to the ratios (fTUj) the statistical uncertainties on fo(q 2 ) become about 5 4- 20%. 
The quality of the results is shown in Fig. El for one of the combinations of quark masses 
used in Ref. ^B]- The points are paired because both K — ► tt and tt — > K transitions are 
considered. 

In order to extrapolate the scalar form factor to q 2 = we have considered three different 
possibilities, namely a polar, a linear and a quadratic fit: 

fo(q 2 ) = /o M °(0)/(l-A ^V), (11) 
h{q 2 ) = rf m >) " (1 + A? m V) , (12) 
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f„(q 2 ) 




o.o 

-0.30 -0.25 -0.20 -0.15 7 0.10 -0.05 0.00 0.05 

(aq) 2 

Figure 3: The form factor fo(q 2 ) obtained from the double ratios HU) for q 2 < g^ ax (full dots) 
and from the double ratio Ffy at q 2 = g^ ax (open dot), for the quark mass combination k s = 
0.13390 and kg = 0.13440. The dot-dashed, dashed and solid lines correspond to the polar, 
linear and quadratic fits given in Eqs. hllMty . respectively. The inset is an enlargement of 
the region around q 2 = 0. 

f (q 2 ) = / M (0).(1 + Ai^ ? 2 + ^ 4 ). (13) 

These fits are shown in Fig. Eland provide values of both /o(0) and the slope Ao, which are 
consistent with each other within the statistical uncertainties. 

Lattice artifacts on /o(0), due to the finiteness of the lattice spacing, start at 0(a 2 ) 
and are proportional to (m s — me) 2 , like the physical SU(3)-breaking effects. Indeed, the 
determination of /o(<?max) * s affected only by discretization errors of 0[a 2 {m s — me) 2 ], because 
the double ratio (J2J) is (9(a)-improved and symmetric with respect to the exchange m s <-» me 
in the weak vertex. In addition, since g^ ax is proportional to (m s — me) 2 , 0(a 2 ) effects in 
the extrapolation from g^ ax to q 2 = also vanish quadratically in (m s — me) . Being in our 
calculation a -1 ~ 2.7 GeV, we expect discretization errors to be sensibly smaller than the 
physical SU(3)-breaking effects. 

The results obtained for /o(0) agree well with a quadratic dependence on a 2 AM 2 , where 
AM 2 = M\ — M 2 , as expected from both physical and lattice artifact contributions. This 
feature is clearly visible in Fig. 01 where the values of /o(0) obtained from the quadratic fit 
()13|). are plotted. 
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Figure 4: Values o//o(0) = /+(0), obtained from the quadratic fit < f73|) . versus (a 2 AM 2 ) 2 . 
The solid line is the result of the linear fit /o(0) = 1 — A (a 2 AM 2 ) 2 where A is a mass- 
independent parameter. 

2.4 Results for the slopes 

Our results for the slope A , extrapolated to the physical meson masses (using a linear 
dependence in the quark masses) and given in units of M 2 + , yield: A{f = 0.0122(22), 
x (un.) = 0.0089(11) and A^ W ° = 0.0115(26). Our "polar" value X^ oL) is consistent with 
the recent determination from KTeV Aq P ° L ' ) = 0.01414 ± 0.00095 [2E| and represents a true 
theoretical prediction, having been obtained before the KTeV result were published. 

As far as the vector form factor is concerned, the values obtained for its slope A + (employ- 
ing the polar parametrization) agree well with the inverse of the squared i^*-meson mass, 
for each combination of the simulated quark masses. A simple linear extrapolation in terms 
of the quark masses to the physical values yields A + = 0.026 ± 0.002 in units of M 2 + , which 
is consistent with the PDG value A + = 0.028 ± 0.002 jB] (obtained from linear fits to the 
vector form factor). The agreement is even better with the recent KTeV result obtained 
using a polar fit: A + = 0.02502 ±0.00037 j2H]- Accurate determinations of the slopes A + and 
A have recently been reported also by ISTRA+ [32] and NA48 (A+ only) (22] • The results 
obtained by means of linear parameterizations are in good agreement with those obtained 
by KTeV. The situation of the quadratic fit to the vector form factor is more controversial. 
For the polar form, which is the one suggested also by the lattice calculation (that explores a 
large range in q 2 ), the NA48 result is in agreement with the one obtained by KTeV, whereas 
this parameterization has not been investigated by ISTRA+. The choice of the polar form 
seems to be the best compromise among the various results, and for this reason it has been 
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used in extracting the experimental values of \V US \ ■ /+(0) presented in Table El 
2.5 Extrapolation of f+(0) to the physical masses 

In order to determine the physical value of /+(0) the lattice results of Fig. H]are extrapolated 
to the physical kaon and pion masses. The problem of the chiral extrapolation is substantially 
simplified if the AG theorem (which also holds in the quenched approximation) is taken into 
account and if the leading (quenched) chiral logs are subtracted. Thus in jTH] the following 
quantity is introduced 

R(M M) ^ + mM K ,M n )-U(0;M K ,M n ) 

R(M K , M n ) (a 2 AM 2 )2 > ( 14 ) 

where /f represents the leading non-local contribution, determined by pseudoscalar meson 
loops within quenched ChPT, calculated in Ref. (TE]. The subtraction of f 2 is a well defined 
procedure being finite and scale- independent. It should be emphasized that this subtraction 
in Eq. (|14|) does not imply necessarily a good convergence of (quenched) ChPT at 0(p 4 ) 
for the meson masses used in the lattice simulation. The aim of this procedure is to access 
directly on the lattice the quantity 1 + /| — /+(0), defined in such a way that its chiral 
expansion starts at 0(p 6 ), independently of the values of the meson masses. After the 
subtraction of /f we expect that, in the range of masses considered in the lattice simulation, 
the ratio R(Mx, M n ) receives large contributions from 0(p 6 ) local operators in the effective 
theory. At the same time, the quadratic dependence on a 2 AM 2 , driven by the AG theorem, 
is already factorized out. Hence R(Mk, M t ) should be a quantity well suited for a smooth 
polynomial extrapolation in the meson masses. It turns out, indeed, that the dependence of 
R(Mk, M n ) on the meson masses is well described by a simple linear fit: 

R^ n -\M K ,M 7T ) = c 11 + c 12 [(aM K ) 2 + (aM 7T ) 2 ] , (15) 

whereas the dependence on AM 2 is found to be negligible. In order to check the stability of 
the results, quadratic and logarithmic fits have been also considered: 

R^ ad -\M K ,M n ) = c 21 + c 22 [(aM K ) 2 + (aM 7T ) 2 } + c 23 [(aM K ) 2 + (aM 7T ) 2 } 2 , (16) 
R«°°-\M K ,M W ) = c 31 + C3 2 log[(aM x ) 2 + (aM^) 2 ]. (17) 

In Fig. it is shown that linear, quadratic and logarithmic functional forms provide equally 
good fits to the lattice data with consistent results also at the physical point. 

Combining our estimate of R(Mk, M n ) at the physical meson masses with the unquenched 
value of f 2 (f 2 = —0.023 |H|), we finally obtain the result [16 j 

ff*~ (0) = 0.960 ± 0.005 stat ± 0.007 syst , (18) 

which has been anticipated in Eq. (j3J). Note that the systematic error does not include an 
estimate of quenched effects beyond 0(p 4 ). Removing this error represents one of the major 
goal of future lattice studies of decays. The error quoted in Eq. ()18j) is mainly due to 
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Figure 5: Comparison among linear < 11,5|) . quadratic and logarithmic \17\ ) fits of the 
ratio R(Mk, M n ) as a function of \a 2 M\ + a 2 M 2 ]. Triangle, square and diamond are the 
values of R(Mk, M n ) extrapolated to the physical meson masses. For illustrative purposes 
we have chosen the case in which a quadratic fit in q 2 is used to extrapolate the scalar form 
factor to q 2 = 0. 

the uncertainties resulting from the functional dependence of the scalar form factor on both 
q 2 and the meson masses. It can be further reduced by using larger lattice volumes (leading 
to smaller lattice momenta) as well as smaller meson masses. We stress again that in our 
estimate of /+(0) discretization effects start at 0(a 2 ) and are also proportional to (m s — 
rrii) 2 , as the physical SU(3)-breaking effects. Thus our result is not affected by the whole 
discretization error on the three-point correlation function, but only by its smaller SU(3)- 
breaking part. Discretization errors on [/+(0) — 1 — f?\ are estimated to be few percent of 
the physical term, i.e. well within the systematic uncertainty quoted in Eq. (j!8|) . For a more 
refined estimate of these effects, calculations at different values of the lattice spacing are 
required. 

3 Semileptonic Yi~ —> uiu decay 

Semileptonic hyperon decays represent the "baryonic way" to a precise determination of 
\V US \. In this Section we present the results of a preliminary lattice study of the vector form 
factor /i(0) for the £~ — > niu decay, based on a strategy similar to the one developed for 
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Figure 6: Results for fo(lmax)' f or hyperon £ — > ntv decays, versus a (M| — M, 



kaon decays. More details on this calculation can be found in Ref. 

The relevant matrix element of the weak vector current can be decomposed in terms of 
the following structures 

<n(p')\urs\E-(p) >=u n (p') | 7 ^i(g 2 ) ~ m^m/ 2 ^ + M n + Mz Mq2) } MP) (19) 
with q = p — p' . As in the case of decays, one introduces the scalar form factor fo(q 2 ): 

fo(q 2 ) = h(q 2 ) + M / M , h(q 2 ) , (20) 

so that /o(0) = /i(0). Note that the SU(3)-symmetric value of /i(0) is given by a Clebsch- 
Gordan coefficient, equal to (—1) for the £~ — ► n transition. 

The value of fo{q 2 ) at g^ ax = (M s — M n ) 2 can be extracted using the double ratio (J7J), 
which now reads as 

^ M) ^^w^m^^) " [/o(gmax)] • (21) 

The results obtained for foiQmax) are shown in Fig. ® where the high precision reached 
(< 0.1%) can be appreciated. 

Given the high accuracy reached at q^ax as well as the closeness of the values of g^ ax to 
q 2 = 0, it is enough to study the g 2 -dependence of /o,i(<? 2 ) with an accuracy of about 10-i-20% 
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(b) versus a 2 q 2 . The dashed and solid lines 
: fit to the lattice data, [see Eq. 



in order to reach the percent precision on /o(0). As in the case of mesons, the standard form 
factor analysis provides values of fi(q 2 ) quite well determined, whereas for fo(q 2 ) one has to 
resort to the double ratios ©, which give access to the quantity fo{q 2 )/fi(q 2 )- In Fig. [7] we 
present the values of fo{q 2 ) and fi{q 2 ) obtained for a specific combination of quark masses. 
The points are paired because both £~ — > n and n — > £~ transitions are considered. 

In order to get the value /o(0) = /i(0), the results for fo(q 2 ) and fi(q 2 ) are fitted by using 
either a monopole or a dipole functional form, 

F (^) tf) = _A_ , F <*. V) = _!_, (22) 

which nicely describe the lattice data, see Fig. In the case of the vector form factor, the 
dipole parameter \J~D agrees with the K* meson mass within ~ 15% accuracy. It should be 
mentioned that, from a theoretical point of view, the dipole parameterization of the form 
factors is only meant to represent, in an effective way, the contribution of several resonances. 
It is found that a dipole behaviour describes well the results of the electroproduction and 
neutrino experiments for AS = transitions. Thus, it might also provide a good description 
of the data in the AS = 1 case (21 Ell- 
in Fig. Owe collect the results obtained for /i(0). They clearly exhibit the nice linear 
dependence expected from the AG theorem. Thus S'?7(3)-breaking effects are resolved with 
a good precision even within the limited statistics used in Ref. |2T] . 

As already pointed out in the Introduction, the evaluation of the leading chiral correction 
f 2 for the baryon form factor /i(0) is more involved with respect to the case of decays. 
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Figure 8: Results for f\(0) versus a 4 (M| — M 2 ) 2 , obtained by performing a monopole (open 
squares) or a dipole (full circles) interpolation in q 2 . The dashed and solid lines are linear 
fits, according to the AG theorem J7|/. 

We believe that the high-precision results obtained for hyperon decays encourage a more 
refined calculation of /2, both in the quenched and unquenched theory, which would permit 
to control the chiral extrapolation of /i(0) at a level of accuracy comparable to the one 
achieved for /+(0) in Ref. [T^J . Work in this direction is in progress. 

4 Conclusions 

We have presented quenched lattice studies of the K — > n and £ — > n vector form factors at 
zero-momentum transfer, /+(0) and /i(0), respectively. Our calculations are the first ones 
obtained by using a non-perturbative method based only on QCD, except for the quenched 
approximation. The main goal is the determination of the SU(3)-breaking effects on /+(0) 
and /i(0), which is necessary to extract |V^ S | from and hyperon decays. In order to reach 
the required level of precision we have employed the double ratio method originally proposed 
in Ref. jHT] for the study of heavy-light form factors. We have found that this approach yields 
a determination of the scalar form factor fo(q 2 ) at q 2 = g^ ax with a statistical uncertainty 
well below 1% for both mesons and baryons. 

A second step is the extrapolation of the scalar form factor to q 2 = 0. This has been 
performed by fitting the accurate results obtained using suitable double ratios of three-point 
correlation functions. The values of /+(0) and /i(0) obtained in this way are determined 
within the percent level of precision, which is the one required for a significant determination 
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In the case of K& decays, the leading chiral artifacts of the quenched approximation, 
represented by /|, have been corrected for by means of an analytic calculation in quenched 
ChPT. After this subtraction, the lattice results are smoothly extrapolated to the physical 
meson masses, and we obtain 

ff*~ (0) = 0.960 ± 0.005 stat ± 0.007 syst , (23) 

where the systematic error does not include an estimate of quenched effects beyond 0(p 4 ). 
Removing this error represents one of the major goal of future lattice studies of decays. 
In this paper we have addressed in some details the impact of the lattice result on the 
determination of \V US \. By considering only the new experimental data on decay widths 
from BNL-E865 jkj, KTeV NA48 [H] and KLOE jZZ], we find \V US \ = 0.2250 ± 0.0021, 
in excellent agreement with CKM unitarity. The results of a preliminary lattice study of 
hyperon decays, based on a similar strategy, have been also presented. 
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